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Abstract 



We construct the nuclear and quark matter equations of state 
at zero temperature in an effective quark theory (the Nambu-Jona- 
Lasinio model), and discuss the phase transition between them. The 
nuclear matter equation of state is based on the quark-diquark descrip- 
tion of the single nucleon, while the quark matter equation of state 
includes the effects of scalar diquark condensation (color superconduc- 
tivity). The effect of diquark condensation on the phase transition is 
discussed in detail. 

PACS: 12.39.Fe, 12.39.Ki, 21.65.+1, 24.85+p, 64.70.Fx, 05.70.-a 
Keywords: Equation of state, Effective quark theories, Diquark con- 
densation. Phase transitions 
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1 Introduction 



The behavior of matter at high baryon density is of great interest in connec- 
tion with neutron stars, the possible existence of quark stars [1], and heavy 
ion colUsions [2] . At normal densities, hadronic matter consists of three-quark 
(nucleon-like [3]) clusters, and the forces between them lead to the familiar 
saturation properties of nuclear matter (NM). At high densities one expects 
a phase transition to quark matter (QM) [4] , but the dynamics of this tran- 
sition, the transition densities, and the equation of state (EOS) in the high 
density regime are largely unknown. In many recent investigations it has 
been found that the QM state at high densities differs from a non-interacting 
quark gas. In particular, there is an instability with respect to diquark con- 
densation, leading to a color superconducting state [5]- [7]. On the other 
hand, diquark correlations are known to play an important role in explaining 
the properties of single nucleons, for example the flavor dependence of the 
structure functions [8, 9]. At high baryon densities one therefore expects a 
phase transition from quark-diquark clusters to a diquark condensed phase. 

Since it is still difficult to obtain unambiguous results for finite density 
matter directly from QCD, effective quark models provide powerful tools to 
investigate the behavior of matter as the baryon density increases. Naturally, 
these models should account for the properties of the free nucleon in terms of 
constituent quarks, as well as for the saturation properties of nuclear matter 
at normal densities. They should also allow for the possibility of a high 
density QM phase, where chiral symmetry is restored and color symmetry 
spontaneously broken. In these respects, the Nambu-Jona-Lasinio (NJL) 
model [10, 11] is an attractive candidate. First, the simplicity of the model 
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allows a direct solution of the relativistic Faddeev equation [12], leading to a 
covariant and successful description of the single nucleon [13] -[16]. Second, it 
has been shown recently [17] that this successful quark-diquark description of 
the single nucleon can be extended to describe the EOS of NM. The essential 
ingredient needed to obtain a saturating EOS was to simulate the effect 
of confinement by introducing an infrared cut-off (in addition to the usual 
ultraviolet one) in order to avoid unphysical thresholds corresponding to the 
decay of hadrons into free quarks [21]. Third, models of the NJL-type are 
able to describe a high density QM state where chiral symmetry is restored 
and color symmetry is broken [5] . It is, of course, desirable to combine these 
three aspects and describe nucleons, NM and superconducting QM in one 
consistent framework. This is the purpose of our present work. 

In this paper we will investigate the role of diquark condensation for the 
phase transition from NM to QM in the NJL model. It is clear that diquark 
condensation will soften the EOS of QM and favor the phase transition. 
The transition densities are therefore expected to decrease with increasing 
pairing strength. We will investigate this relation, and discuss a possible 
scenario for the phase transition. In the course of the calculation we will 
face the problem of whether or not those parameters of the model which 
have been determined exclusively by the properties of single nucleons and 
NM, like the value of the infrared cut-off, the strength of the qq interaction 
in the vector meson channel, and the strength of the qq interaction in the 
scalar diquark channel, should be taken over to QM. In order to discuss these 
points, in particular the strength of the qq interaction in the vector meson 
channel, we will also consider the uj meson mass in the medium. Another 
purpose of this work is to present a consistent formulation which accounts 
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for the three aspects discussed above - i.e., the single nucleon, normal NM 
and superconducting QM. 

In sect. 2 we will first explain the model in the framework of the Nambu- 
Gorkov formalism [22, 23], and then use path integral methods to describe the 
single nucleon as well as the EOS for NM and QM. In sect. 3 we will present 
present numerical results, and discuss the role of diquark condensation in the 
phase transition. Finally, we summarize our results in sect. 4. 

2 NJL model in the Nambu-Gorkov formal- 
ism 

The Lagrangian density of the flavor SU(2) NJL model has the general form 
C — ^{i^ — m) i> + J2a (i^^a'tf) , where m is the current quark mass, and 
the chirally symmetric 4-fermi interaction is characterized by the matrices Fq, 
in Dirac, flavor and color space, and the corresponding coupling constants 
Ga- Hadronization techniques to describe the single nucleon and NM, as 
well as the description of superconducting QM, are simplified very much if 
one makes use of the Nambu-Gorkov formalism [22, 23]. In this method, the 
Lagragian density is re-written identically in terms of the fields ^ 

where C = ^70 72, as follows: 

C^^(i^-m)^ + J2Ga (^F^^M*)' . (2.2) 

aM 

"'^ Since there is freedom under unitary transformations to define the field the presence 
of the matrices C and T2 in the lower components is a matter of definition. We include 
them here because they lead to simplifications in some of the following formulae. 
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Wc denote the matrices in the Nambu-Gorkov space by o"m = {<^o, o";?) for the 
interactions in mesonic {qq) channels, and by an = (o"i,cr2) for the diquark 
(qq) channels. Here ctq is the 2x2 unit matrix, and (7i, (J2, (73 are the usual 
Pauli matrices. It is easy to verify that the only non-zero terms in (2.2) are 
those which satisfy 

CT2a2AlC-'a2r2 = Aa, (2.3) 

where = FqCTm- Therefore, for a given Fq,, either M — or M — 3 
contributes in (2.2). 

Any Lagragian density of the form (2.2) can further be re-written into a 
Fierz symmetric form [12], where the interactions in the diquark channels, 
characterized by ao, also appear. For these terms to be non-zero, the same 
condition (2.3) holds with A^ = TaCTn- The coupling constants appearing in 
the Fierz symmetric Lagrangian density for the various meson and diquark 
channels are linear combinations of the G^s appearing in the Lagrangian 
density (2.2). For the purposes of this work, we will need only the following 
terms: 

- (^7^(73^)' + G',(^^75a^/3„^)', (2.4) 

which are the interactions in the ( J^, T) = (0+, 0), (0~, 1) and (1~, 0) color 
singlet meson channels, and the (0+, 0) color 3 diquark (scalar diquark) chan- 
nel. The color 3 matrices are defined in terms of the antisymmetric Gell- 
Mann matrices by Pa = y^3/2 A^, where a = 1, 2, 3 corresponds to A = 2, 5, 7. 
A sum over D — 1,2 and a = 1, 2, 3 is imphed in the last term of Eq.(2.4). 
The 3 coupling constants in Eq.(2.4) will be treated as free parameters. 
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The values of the ratios 

rs = -pr^ rcu = -pT- (2-5) 

reflect the form of the original interaction Lagrangian density (2.2). 

The 4-fermi interactions, Gq, (^^A^^^ , in Eq.(2.4) can be eliminated in 
favor of bosonic auxiliary fields ba as usual ^ by subtracting terms of the form 

(2Ga^Aa^ + /4:Ga, which are formally equal to zero if we impose the 
constraints (equations of motion) dC/dba = 0. In this way we obtain the 
partially bosonized Lagranian density in the form: 

jC = ^' (i^ - m - S - i75T ■ 7ro-3 - ^0-3 - i^^ao^Dafia) ^ 

where the real bosonic auxiliary fields are defined by 

E = -2(7^**; TT = -2G'^*i75T(T3*, (2.7) 

= 2G,^ 7^(73*; ^Da = -2Gs^ti^(3aaD-^. (2.8) 

In the two terms of (2.6) involving the diquark fields, a sum over D = 1, 2 
and a = 1, 2, 3 is implied. 

3 Nucleons, nuclear matter, and color super- 
conducting quark matter 

In this section we will use the Lagrangian density (2.6) to construct the single 
nucleon as a quark-scalar diquark state. We then use it to obtain the EOS 
of both NM and QM in the Hartree approximation. 

^This procedure can also be formulated in the path integral formalism, as will be done 
in sect. 3 for the nucleon auxiliary field - see Eq.(3.9). 
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3.1 Nucleons and nuclear matter 

The EOS of NM, based on the quark-diquark description of the single nu- 
cleon, has been derived in Ref.[17] by using a hybrid model to evaluate the 
expectation value of the quark Hamiltonian in the NM ground state, similar 
in spirit to the model of Guichon and collaborators [3, 24]. Since it is one 
of the aims of this paper to present a consistent formulation of the single 
nucleon, NM and QM in one framework, we will follow the path integral 
formalism based on the quark Lagrangian density (2.6). 

The general method to obtain the NM EOS is as follows: One first uses 
hadronization methods [25, 14] to derive an effective hadronic Lagrangian 
density of the cr-model type from the quark Lagrangian density, and intro- 
duces a chemical potential for the nucleons. Then one can use any of the fa- 
miliar approximation schemes for relativistic many- nucleon systems [26, 27]. 
In this paper we will limit ourselves to the Hartree approximation. For this 
purpose it is sufficient to derive the nucleon kinetic energy term in the pres- 
ence of the mean meson fields, the c-number terms associated with the mean 
fields, and the "trace-log term" which emerges from the integration over the 
quark fields. Therefore, although the method outlined below can be used 
to obtain more complete expressions for the effective hadronic Lagrangian 
density, we will concentrate here on the nucleon kinetic energy term in the 
presence of mean scalar and vector fields. 

We introduce the (color singlet) nucleon fields and sources by adding the 
following terms to the Lagrangian density (2.6): 



$ + $ 



+ ^ i a: r ^ 



(3.9) 



Here the complex diquark fields are defined by A^ = Ai^ — i^2a and A* 



Ala + i^2a, and the Nambu-Gorkov fields, A'' and are expressed in terms 
of the nucleon field, N ^ and the nucleon source, 0, in a similar manner to 
(2.1) for the quark fields. In the generating functional, the integrations over 
$ and $ give functional 5-functions which show that the nucleon auxiliary 
fields must satisfy the constraints N — —t/ja^a/V^ and N — — ■0„A*/-\/3. 
The integrations over the nucleon fields then give unity - i.e., the terms (3.9) 
do not change the generating functional of the theory. 

After integration over the quark fields (see Appendix A) we obtain the 
effective action ^ 

= --Trln5 +Jdx[-^ AG^r^^J ^'"'^^ 

Here 

S-^ ^j^-m-T,- i75T • 77(73 - 0(^3 - «75 1^ ^q'^" j , (3.12) 

and S in (3.11) is the inverse of this. The meson and diquark fields are now 
the dynamical variables, while the nucleon field is still an auxiliary variable. 

In NM the scalar and vector meson fields have non-zero expectation val- 
ues, and accordingly one shifts the meson fields to separate these c-number 
parts from the fluctuation parts over which one has to integrate. Since we will 
describe the nucleon as a quark-diquark state (omitting mesonic fluctuations 
for the present), and the NM EOS in the Hartree approximation, we omit 

^In the action, the symbol Tr inchidcs also a functional trace. In order to avoid too 
many symbols, fields (propagators) will be frequently considered as vectors (matrices) in 
function space, besides their structures in color, flavor, and Dirac space. The notation <Seff 
will be used generically for the effective action, irrespective of the space of fields where it 
operates. 
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the fluctuation parts of the bosonic flelds from now, except for the scalar 
diquark field which will be used to construct the nucleon. That is, from now 
E and o;^ will denote only the classical parts of the scalar and vector fields, 
which are expressed in terms of the quark fields ip, ip as ^ 

E = -2G'^(NM|V^V|NM) = M-m, = 2G^{NM\i^ Yi^lNM) (3.13) 

where |NM) is the NM ground state. Since the pion field is assumed to have 
zero expectation value in NM, it will be omitted from now on. 

The integration over the diquark fields can no longer be done exactly. 
Here we resort to the stationary phase approximation, where only the terms 
quadratic in the diquark fields are retained in the quark determinant - i.e., 
n-point diquark interaction terms (n=4,6,8. . .) are omitted. After performing 
the Nambu-Gorkov trace in (3.10), this amounts to the replacement 

-'-TrlnS-'- Jd'x^^-itrlnS^' + A^D^'A:. (3.14) 

Here we introduced the following quantities, which we write down in momen- 
tum space for further use: 

n,(*) = 6i I J^tiB [7550(5)7550(9 - *:)] (3.16) 



""^The quark mass parameter M is considered as a variable in the effective potential, and 
its physical value, which is the solution of the gap equation at density p > (p = 0), is 
denoted as M* (Mq). The same notation will be used for the nucleon mass Mjv = MAr(M) 
and the gap A, i.e., their physical values are = Mn{M*) and A* at finite density, 
and Mjvo = Mjv(-Mo) and Aq = at zero density. 
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Since we wish to construct the nucleon as a quark-diquark state, we will 
retain the diquark fields to all orders in the first term of (3.11) in order to 
derive the nucleon propagator in the ladder approximation to the Faddeev 
equation. To evaluate this term, we need the form of S, which is easily de- 
rived by inverting (3.12) - see Appendix A. Since only terms with the same 
number of A's and A*'s survive after the diquark integration, because of 
baryon number conservation, the non-diagonal parts of S in Nambu-Gorkov 
space, which involve unpaired fields A or A*, do not contribute. After per- 
forming the Nambu-Gorkov trace, we therefore find that the first term in 
(3.11) becomes effectively (see Appendix A) 



After the integration over the diquark fields the nucleon field becomes a 
dynamical variable, and this term gives the Faddeev propagator. In order to 
simplify the formulae, we will refer here only to the "static approximation" to 
the Faddeev equation [28] , which will be used for the numerical calculations 
in this paper. It is obtained by replacing Sq, which appears between the 
diquark fields in the denominator of (3.18) and describes the quark exchange 
between the diquark and the spectator quark, according to 5*0 —1/M. ^ 

The integration over the diquark fields is performed as usual by treating 
the sum of the second term in (3.14) and the term in (3.18) proportional to 
AA* as the free diquark action 




(3.18) 




(3.19) 



^Only in this approximation is the nucleon auxihary field a local one. 
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and the remaining part of (3.18) as the interaction part 




^ {AM (KPa) 



SoA* <p. 



(3.20) 



(A„,/3„0 (A*/3„) 



After introducing an auxihary source term A*Ja + Ja^a into the Lagrangian 
density to perform the diquark integration, we obtain, with Sd = Sdo + S^j: 



Expanding this expression in powers of the nucleon sources gives the nucleon 
propagator as the coefficient of 00, and interactions between miclcons as the 
higher order terms ^. Here we concentrate on the nucleon kinetic term. Since 
the interaction part, Eq.(3.20), contains all powers of the diquark fields, many 
contractions emerge when the derivatives indicated in (3.22) are carried out, 
even if we restrict ourselves to the nucleon kinetic term ex 00. A subset 
of these contractions constitutes the ladder approximation to the Faddeev 
equation (see Appendix A), which leads to 



where the nucleon (quark-diquark) propagator in the ladder approxima- 
tion is expressed in terms of the quark-diquark t- matrix (Tat) in momentum 

^An interesting example is the second exponential factor in (3.21), which describes 
interactions between nucleons proceeding via quark exchange processes. 




(3.21) 



(3.22) 




(3.23) 
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space by 



(3.24) 



(3.25) 



Tn{p) 



-3 1 



(3.26) 



M 1 + An^(p) 



The quark and diquark propagators, -Do and So, appearing in the quark- 
diquark bubble graph, (3.25), have been given in Eqs.(3.15) and (3.17). ^ 
As explained in sect. 2. 1.1 of Ref.[17], Tjv has a pole at = M^, where 
Pn — P^ ~ 3a;''. If the total system under consideration (nuclear matter) is at 
rest, the space components of the vector field vanish, and the nucleon positive 



energy spectrum becomes eiv(p) = En{p) + 3a;°, with En(p) — Jp^ + M%. 



Inserting these results obtained above into the effective action of Eq.(3.10), 
(3.11), and finally integrating over the nucleon sources (see Appendix A), we 
obtain the effective action for the nucleons in the mean field (Hartree) ap- 
proximation: 



The three "trace-log terms" in (3.27) are the results of the integrations over 

the quark fields, the diquark fields, and the nucleon sources, respectively. 

The first term in (3.28) is the kinetic term of the nucleon, and the rest are 

the c-number terms associated with the mean scalar and vector fields. 

^The diquark propagator, Dq, is related to the quantity Ts of Rcf.[17] by Dq = its- The 
t-matrix (3.26) is related to T{p) in Eq.(2.7) of Ref.[17] by Tn{p) = -T{p). We also note 
that the factor 2Gu in Eq.(2.8) of the present paper was not included in the definition of 
the field w'^ in Ref.[17]. 





(3.27) 



(3.28) 
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We now introduce a chemical potential (n) for the nucleon into the kinetic 
term of (3.28). Before that, however, we note that if one had continued 
and integrated exp (iSes) over the nucleon fields as well, the last term in 

(3.27) would cancel with the result of the integration, and one would get an 
expression which could be obtained directly from (2.6) without any reference 
to nucleons. This is obvious, since we introduced the nucleon auxiliary fields 
in (3.9) in order that the generating functional would be unchanged. This 
observation is important, however, because it means that, as a consequence of 
the composite nature of the nucleon, there are no nucleon vacuum loop terms. 
That is, there are no nucleon loop contributions to the effective potential 
unless we introduce a chemical potential into the kinetic energy term of 

(3.28) . The vacuum loop terms are exhausted already by the quark and 
diquark loops in (3.27), which describe the polarization of the quark Dirac 
sea and the effects of quark-quark correlations (so called, ring contributions) 
in the quark Dirac sea, respectively. 

Concerning the ring contributions (corresponding to the second term in 
Eq.(3.27)), we note that meson loops, which were omitted from the outset 
in our simplified treatment, would give very similar terms. A consistent 
treatment of the ring contributions would therefore require both meson and 
diquark loops. This is beyond the scope of the present approximation, where 
we have retained the diquark fields solely in order to construct the single 
nucleon as a quark-diquark state. In the framework of the Hartree approxi- 
mation for the EOS, we therefore drop the second term in Eq.(3.27). 

One can rewrite the nucleon kinetic term by introducing a renormalized 
nucleon field (N), such that, near the pole, the renormalized propagator (Gn) 
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behaves as G 



N 



M 



N- 



The pole part of the nucleon kinetic term 



including the chemical potential then becomes 



(3.29) 



where fj,* = fj, — 3 uj^ . 

To continue the derivation of the effective potential {V) for NM in the 
path integral formalism, we finally have to integrate over the nucleon field and 
express the result as Seff — — J d'^xV, where is a function of the classical 
fields M and Combining the result of the integration of the kinetic term 
(3.29) with the last term of Eq.(3.27), we get the nucleonic contribution to 
the effective potential^ 



VN = t I T^tr In - + //*70) - (/x = term) 



(27r) 



- 2><2 



= 



-4 



d^p 
(2^ 

d'V 
(2^ 



In 



tro In (-^ -Mn + //*7°) - + //*7°) 



(3.30) 
(// = term) 



+ In 



En{p? 



Po 



ENipf 



-3a;°) 
(3.31) 
(3.32) 



which is nothing but the familiar contribution to the effective potential arising 
from the nucleon Fermi motion. Adding the quark loop contribution (the 
first term in (3.27)) and the mean field terms of (3.28), we finally obtain the 

iZf^N and Gjv = ZnGn with Zn = 



^Mathematically, wc write N = 

-(M/3) [(5njv(pAr)/a^jv)^„=M«]~'- 

^For the explicit evaluation of the Dirac determinant, see Appendix B for the case 
A = 0. The trace indicated by the symbol tr does not include the functional trace. The 
subtraction of the = term in (3.30) arises from the last term in (3.27), and the factor 
2 in the second line comes from the isospin trace. 
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effective potential for NM in tlie form 




(3.33) 



where is given by (3.32), and the vacuum term is 



J (2^ A;2 - Mi 



+ 



(M - mf _ (Mo - mf 



(3.34) 



Here we subtracted the zero density contribution (M = Mq) . The values of 
loq and M for fixed // are determined by the equations^" dV / duj=dV / dM — 0. 
The pressure as a function of jji is then obtained as P = —V. The baryon 
density follows from p = —dV/dp, and the energy density is given by S{p) = 
V + pp. 

The effective potential (3.33) differs from the familiar expression in the 
linear sigma model for point-nucleons [29] only through the dependence of the 
nucleon mass Mjv on the mean scalar field (or, equivalently, on M). Since the 
scalar field now couples to the quarks in the nucleon instead of an elementary 
nucleon, the function M^{M), which is determined by the pole position of 
the quark-diquark t-matrix (3.26), is in general a non-linear function of M. 
If it has a large positive curvature (scalar polarizability), the binding energy 
per nucleon saturates, as has been discussed in detail in Ref.[17]. 

3.2 Color superconducting quark matter 

To obtain the EOS of QM in the Hartree approximation, we introduce a 
chemical potential for quark number pg = p/3 into the Lagrangian density 

^°In principle, the quantity which has to be minimized is the effective potential after 
Wick rotation. In this case the last term in Eq.(3.33) changes its sign, and the physical 
value of loq is a minimum of V. 
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(2.6), and assume that the fields S, u;° and An have finite expectation val- 
ues. The choice I? = 1 for the diquark field Ad^ corresponds to a special 
choice of the phase and breaks the C/(l) symmetry, while the choice a — 1 
breaks the color SU{3) symmetry down to SU{2) Since we do not con- 
sider fluctuations of the bosonic flelds in this work, it is sufficient to use the 
following Lagrangian density: 

Here 

S-^^i^-M + 7° (73 - i75 (Ti A (3.36) 

where /j,* — Hg — Uq, and the meson mean fields are defined in terms of the 
ordinary quark fields ip, ip by (compare with Eqs. (2.7), (2.8), and (3.13)) 

M = m-2G^< QMIV^V'IQM >, uo = 2G^(QM|V^V|QM), 

A = -2G,^(QM|(v^i75/3iCT2V^^-V'^C-V2i75/3iV')|QM). (3.37) 

The poles of the quark propagator, S, are easily be determined by invert- 
ing Eq.(3.36). There are 4 poles, at po — ±-sJ{Eq(p) ± //*)2 -|- A^, each with 
degeneracy 8, and 4 poles at po — ±\Eq{p) ± each with degeneracy 4. 



Here Eg{p) = y/pH^JP, and A = 




^"'^Thcrc arc 5 different Goldstone bosons corresponding to these broken symmetries. 
If we start from the (ai /3i) mode (An), which characterizes the ground state, the [/(I) 
phase rotation mixes in the ((72 /3i) mode, and the color SU{3) rotation mixes in the (cri ^2), 
{o'l /^a) and ((T2 Pa) (a = 1, 2, 3) modes, which correspond to the Goldstone bosons. Note 
that the (cr2 (3i) mode is mixed into the ground state by both phase and color rotations. 
The quark-quark t-matrices in these channels have zero energy poles. For comparison, we 
note that the familiar Goldstone pion in the present Nambu-Gorkov formalism emerges as 
a cr3-mode in the (0~,T = 1) color singlet channel, which is mixed into the cri-mode in 
the (0+,T = 0) color singlet channel of the ground state by the chiral rotation. 
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Integration over the quark fields gives a contribution — -Trln^*"^ to the 
effective action, which can be evaluated (Appendix B) with the result 

-^Trln^-^ ^- Jd'^x 

""^^ / " ^^^^^^ + - A^) + 2 In {pI - {E,{p) - - A') 

+ln {pI - {E, + f,;)') + In {pI - {E, - f.^)] . (3.38) 

If we separate the contributions surviving for A = 0, we finally obtain the 
following form of the effective potential for QM (compare to (3.33) for the 
NM case): 

V^V^ + Vq + Va-^ (3.39) 

where Kac has the same form as in nuclear matter (see Eq.(3.34)), Vq is the 
contribution of the quark Fermi motion given by 

Vq = -12 / - E,{p)) - E,{p)) , (3.40) 

and Va is the contribution arising from the finite gap: 



Va = 4t 



pI - {E,{p) + ii*Y pI - {Eg{p) - II*) 



A^ 

The values of a;o, M and A for fixed chemical potential are obtained by 

solving the equations dV/dujQ — dV/dM — dV/dA — 0. (Concerning the 

dependence on cuq, see the last footnote in subsect.3.1.) The pressure, baryon 

dV 

density and energy density are then obtained from P — —V, p — — and 

a/i 

£ = V + lip, where p = 3pq is the chemical potential for baryon number as 
before. 
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4 Numerical results 



In the numerical calculations discussed below, the proper time regularization 
scheme will be used. In this scheme, one performs a Wick rotation of the 
energy variable of the loop integral, and then introduces the replacements 



where A depends on the momenta and, possibly, Feynman parameters, and 
Auv and Air denote the UV and IR cut-offs, respectively. 

As discussed in detail in Ref.[17], the vahdity of the "static approxima- 
tion" to the Faddeev equation, which leads to the contact-type quark-diquark 
interaction oc 1/M in the t-matrix (3.26), breaks down as M decreases with 
increasing density, and leads to large deviations from the exact Faddeev re- 
sults for the function Mjv(M). On the other hand, it was shown that if one 
fixes the quark-diquark interaction to its value at zero density (1/M l/Mo 
in the denominator of Eq.(3.26)), the exact Faddeev results are qualitatively 
reproduced. Furthermore, if one uses an interpolating form of the quark- 
diquark interaction (1/M (l/Mo)(Mo + c)/(M + c), with c = 0.7 GeV, 
in the denominator of (3.26)), the agreement is very good. We will use this 
interpolating form in the following calculations. 



In order to provide a basis for our later discussions on the NM^ QM phase 
transition, we first reproduce in Figs. 1 and 2 the results obtained in Ref.[17] 
for the nucleon mass as a function of the scalar potential, M — Mq, as well as 




(n > 1) 



(4.42) 



4.1 EOSofNM 
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the binding energy per nucleon as a function of density. The parameters used 

in this calculation are shown in Table 1. For the choice Ajr — there is an 

unphysical threshold for the decay of the nucleon into a quark and a diquark, 

and the function Mn{M) cannot develop a sufficiently large curvature (scalar 

polarizability) , which is required to stabilize the system. In contrast, for the 

case Air = 0.2 GeV the nucleon mass exceeds the would-be threshold for 

large scalar potentials, leading to a sufficiently large scalar polarizability and 

to saturation of the binding energy The effective masses of the nucleon 

and the quark for the case Air — 0.2 GeV are shown by the solid lines in Fig. 

3. We observe that there is only a limited tendency toward chiral restoration 

in NM. The pressure of NM will be discussed later in connection with the 

phase transition to QM. 

Let us discuss here another aspect which provides a further illustration of 

the importance of avoiding unphysical decay thresholds, namely the mass of 

the a;-meson in NM. We take the simplest case of an cu- meson at rest (q = 0), 

and consider the pole position qq = M* of the effective NN interaction in the 

t-channel with quantum numbers 1+. The effective NN interaction for the 

case q — qQ — Q\s the familiar Landau- Migdal interaction, which was derived 

in Ref.[17] directly from the expression for the NM energy density. The part 

associated with the spatial components of the exchanged a;-meson has the 

^^Wc recall from Ref.[17] that the only free parameter for the calculation of the NM 
EOS, rui, has been adjusted so that the binding energy curve passes through the empirical 
saturation point. That is, with the limited number of parameters in this simple NJL 
model, it is not possible to ensure that the calculated saturation point agrees with the 
empirical one. 



20 



form 



p'p 



(4.43) 



where Ep = ypp + M^. The part involving ISG^j^ is a density dependent 
N-loop contribution - i.e., the Fermi average over the Z-graph for external 
vector fields with zero momentum. (This is the same type of medium correc- 
tion which reduces the enhanced isoscalar magnetic moments in relativistic 
many-nucleon theories [30].) 

Since the nucleon remains a rather massive object even at finite density 
(see Fig. 3), we can expect that for the case of finite energy transfer this 
N-loop contribution will not depend very strongly on Qq, and wc therefore 
approximate it here by its value at go = 0. Then, for finite qo, the only 
additional contribution arises from the exchange of qq pairs in the t-channel, 
described by the transverse part of the polarization (bubble graph) r!;^ - see 
Appendix C. The constant in the denominator of Eq.(4.43) is then replaced 
according to 



The two terms on the r.h.s. of Eq.(4.44) work in opposite directions for 
increasing density; the N-loop term represents the effect of the Pauh exclusion 
principle acting on the nucleons and tends to make the 6i;-meson heavier, while 
the qq excitation piece tends to make the a;-meson lighter (since the mass of 
the qq pair is reduced). 

Using the density dependent quark mass in NM, as shown in Fig. 3, and 
the value of r^^ from Table 1, we obtain the result for M* shown in Fig. 4. 

^"^Scc Eq. (2.50) of Ref. [17] . We note that the factor i appearing there in the denominator 
is a misprint. The correct expression is (4.43). 




lSG^-^ + 2G^U^{qo). 



'F 



(4.44) 
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First we note that has been adjusted in Ref.[17] so that the binding energy 
curve passes through the empirical saturation point. From Fig. 4 we see that 
this leads to an u; meson mass of 827 MeV at zero density, which is larger than 
the experimental value of 783 MeV, but nevertheless reasonable. Second, the 
CO meson mass decreases by about 80 MeV at normal NM density and by 
about 120 MeV for p ~ 0.5 fm~^. For higher densities it increases slightly 
because of the increasing importance of Pauli blocking In comparison 
with our earlier NJL model calculations, we note that the overall picture for 
the a;-meson mass is similar to that for the cr- meson mass (see Fig. 11 of 
Ref.[17]), and also to the nucleon and quark masses shown in Fig. 3. That 
is, the density dependence is rather mild, and there are no sudden changes of 
the behavior at high density. Third, we note that the cu meson mass shown 
in Fig. 4 is always above the would-be qq threshold. If we set Air = 0, 
leaving the other parameters unchanged, wc would still get a bound state 
at zero density (M(^ = 780 MeV), but this state would already dissolve at 
p ~ 0.04 fm"", because of the decreasing qq threshold. Once again this shows 
the importance of avoiding the unphysical threshold in the NM calculation. 

4.2 EOS of QM without the effect of diquark conden- 
sation 

The effective quark mass in QM is shown as a function of the density by the 
dashed line in Fig. 3. The parameters used in this calculation are the same 
as for NM, see Table 1. (We note that the value of r^j has no influence on 

^"^In view of the interest in searching for possible bound (x;-nucleus states [18]. it is 
interesting to note tliat these results are very close to those found in the quark meson 
coupling model [19], with the w-meson being somewhat less bound than it would be in 
Walecka-type models [20]. 
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M*.) The sharp decrease of M* in QM around 0.3 fm~'^ reflects the well 
known chiral phase transition which would happen for the case m — 0^^. 
We see that the NM and QM cases are quahtatively very different and one 
should not treat ordinary finite density matter as QM. 

Fig. 5 shows the variation of the pressure as a function of the chemical 
potential for NM (solid line), and for QM with the same model parameters 
as for NM (dashed line labeled by = 0.37), as well as for the choice r^^ = 
(dashed line labeled by = 0). The case of zero density corresponds to the 
point (P, /x) = (0, Mm) for NM, and to (0, 3Mo) for QM, where M^o = 940 
MeV and Mq = 400 MeV are the zero density values of the nucleon and con- 
stituent quark masses. Starting from this point, the density increases along 
the lines. The vacuum solution P — 0, which corresponds to a minimum 
of the effective potential for fi < M^q in NM, and for fj, < 3Mq in QM, is 
not shown in this plot. For a qualitative understanding of our following dis- 
cussions it is helpful to recall that attractive effects, which soften the EOS, 
make the P(/x) curves steeper, and repulsive effects make them flatter. 

Let us first summarize some well known points concerning the phase tran- 
sitions which occur within the NM and QM phases separately. In the region 
below the saturation density, the NM EOS shows the familiar behavior of 
a first order gas-liquid phase transition, which is shown in more detail in 
the insert of Fig. 5. The pressure decreases with increasing density until 
it reaches a minimum at some density p^, and then it increases and passes 
through P = at the saturation density po- In other words, for densities 
between pc and po there are 3 extrema of the effective potential as a function 

Since in the numerical calculations of this paper we always adjust the parameters 
so as to reproduce the experimental pion mass, we refer to this kind of behavior as the 
"would-be chiral phase transition" , or simply the "chiral phase transition" . 
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of M, corresponding to (i) a maximum at negative pressure corresponding 
to unstable NM droplets, (ii) the NM phase, which has negative pressure for 
p < Po (quasi-stable NM droplets), and positive pressure for p > po (stable 
NM), and (iii) the vacuum phase {P — 0), which is the ground state (state 
with the largest pressure) for p < po. 

The EOS of QM for the case = has a bound state, and in this case 
the behavior is qualitatively similar to that of NM, discussed above. The 
gas-liquid phase transition of NM now corresponds to the familiar first order 
chiral phase transition in QM, where the unstable branch with decreasing 
pressure is characterized by massive constituent quarks (broken chiral sym- 
metry), and the metastable {P < 0) and stable {P > 0) branches with 
increasing pressure correspond to the chirally restored phase. The first order 
phase transition occurs at P = in this case, but if one increases r^^ the curve 
will cross itself at finite P, and a further increase of leads to a chiral phase 
transition of second order, as shown for r^^ = 0.37 in Fig. 5. In this case there 
is only one minimum of the effective potential as a function of M for any 
given p,, and at some p the maximum at M ~ turns into a minimum. The 
density increases along the dashed line for — 0.37, starting at the point 
where // = 3 Mq = 1200 MeV and P = 0, and at // = 1500 MeV (p = 0.41 
fm~^) the quark effective mass is already as small as M* = 53 MeV. 

The branch of the QM EOS corresponding to (almost) massless quarks is 
essentially the EOS of a free massless quark gas [4]. This is demonstrated by 
the dotted lines in Fig. 5, which show the EOS for a massless quark gas for two 
choices of the bag constant B: -Bnjl = 139.7 MeVfm"^ is the height of the 
Mexican hat vacuum potential (3.34) at M = 0, and -Bmit = 57.5 MeVfm"^ 
is the value used in the MIT bag model. By comparing the dashed line for 
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r^o = with the dotted hne for B = -Bnjl, we see that, in the region which is 
relevant for a discussion of the NM^ QM phase transition, one could neglect 
the spontaneous breaking of chiral symmetry from the outset and start with 
massless quarks. 

By comparing the solid hne with the dashed line for the case Vu, — 0.37 
in Fig. 5, we see that there is no NM — > QM phase transition if we use the 
same model parameters to describe the two phases. That is, NM would be 
the ground state for all densities. It will be demonstrated later that this 
unsatisfactory situation cannot be improved by including the possibility of 
diquark condensation. (A stiffer EOS for NM would certainly come closer 
to the QM result at high densities, but one never gets a crossing for any 
reasonable NM EOS.) On the other hand, the dashed curve with — 
is much closer to the NM curve, although still there is no phase transition. 
However, it will be shown below that this situation changes immediately if 
one allows for the possibility of diquark condensation. We conclude from 
Fig. 5 that it is not possible to get a phase transition from NM to normal 
QM by using the same parameters for the description of the two phases, in 
particular the same strength of the vector interaction. This point was also 
noted in a recent investigation [32]. However, if one uses r^; = in QM 
and further considers effects which soften the QM EOS, like diquark or pion 
condensation, there could be a phase transition. 

Based on the preceding discussion, we are led to consider the question of 
whether it is reasonable to use the same model parameters for NM and QM 
or not. Two among the parameters hsted in the second column of Table 1 

^^We also note that the dotted hne labeled hy B = Bmit in Fig. 5, which has only little 
to do with the NJL model, crosses the NM line in the high density region, but there is 
also an unphysical crossing in the low density region. 
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have been fixed by considering only the properties of NM, namely Air and r^^. 
First, one needs a finite IR cut-off in order to describe a stable NM state, but 
the results are rather insensitive to its actual value as long as Air > 0.1 GeV. 
The introduction of a finite Air takes into account one aspect of confinement 
physics, which is important for single hadrons and NM. However, there is 
no reason a priori to introduce an IR cut-off in QM. On the contrary, one 
might expect that quark decay processes are welcome in this case, and that 
the choice Air = would be more reasonable. Second, one also needs a finite 
value of in NM in order to reproduce the saturation properties. Once 
again, there is no compelling reason to use the same value in QM, and from 
Fig. 5 we see that the QM EOS depends very sensitively on the actual value 
of r^. 

In order to shed some light on the question of which values of Air and 
arc reasonable in QM, we return to the problem of the uj meson mass as an 
example, and investigate whether or not there is a pole corresponding to the 
u! meson in the qq t-matrix in QM. Contrary to the NM case discussed in the 
previous subsection, we now have a chemical potential for the quarks, which 
leads to a modification of the qq bubble graph arising from Pauli blocking. 
(The detailed formulae are given in Appendix C.) Fig. 6 shows the situation 
for the choice Air = 0.2 GeV in QM. The would-be threshold, which is 
'^\Jp'f + M*"^^ is shown as a function of density by the dashed line, and the 
pole positions for q = are shown for various values of r^^ by the solid lines 

We see that poles exist even for very small values of the coupling constant. 

^^Thc density dependence shown by the dashed hne reflects the behavior of the chemical 
potential /x*, see the dashed line for r,^ = in Fig. 5. The decrease of M* with increasing 
density lowers the threshold, while the effect of Pauli blocking raises the threshold. In 
the low density region the first effect dominates, while in the high density region, where 
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Wc expect that the existence of a hadron pole in QM for all densities is not 
restricted to the case of the lo meson shown here, but can be ascribed to the 
finite IR cut-off. This observation indicates that it is not reasonable to use 
a finite IR cut-off in QM. (Another argument in favor of the choice Air = 
will be given in the next subsection.) 

Fig. 7 shows the situation for the case Air = in QM. We see that for 
the choice = 0.37 there is again a pole for all densities, which now is below 
the threshold. If the value of is decreased, the pole disappears in the low 
density region, but at high densities it re-appears because of the increasing 
threshold. For r^; = 0.17 the pole exists only for p > 0.8 fm~^. Since the 
appearance of a meson pole in high density QM, which was not present in the 
vacuum, is physically unreasonable. Fig. 7 indicates that one should really 
set = in QM. 

Based on these discussions, we will use Air = and r^^ = in QM in 
our following discussions. The corresponding parameter set is shown in the 
third column of Table 1. As compared to the case Air = 0.2 GeV, there are 
small changes in m, and Auv in order to reproduce the same input values 

= 140 MeV, = 93 MeV and Mq = 400 MeV. We will investigate the 

dependence of the QM results on the pairing strength, and therefore we treat 

Ts as a free parameter in QM. (Some numerical results obtained by using the 

same values of Air and/or r^^ as for NM are shown in Appendix D.) 

the quarks are practically massless, the Pauli effect leads to a continuous increase of the 
threshold. Figs. 6 and 7 show that the density dependence of the pole positions follows 
the behavior of the dashes lines 
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4.3 EOS of QM including the effect of diquark conden- 
sation 

The effect of diquark condensation on the QM EOS is represented by the 

piece (3.41) of the effective potential. In order to discuss the nature of the 

phase transition to the color symmetry broken phase, it is useful to look at 
dV 

the gap equation — — = for fixed M. Besides the trivial solution A = 0, 
the nontrivial solution is obtained by solving the equation 

1 = /(A'), (4.45) 

where the explicit form of the function /(A^) in the proper time regulariza- 
tion scheme is 

12G 



f(A')^^ dko edk{F^ + F^), (4.46) 

TT"^ JO JO 



where 



^0 + 4 

with = {Eq{k) ± /x*^ + A^. It is easy to show that /(A^) is a monoton- 
ically decreasing function of A^. For the case A/r > 0, its value at A^ = 
is finite, and therefore we have one unique nontrivial solution if /(O) > 1. 
In other words, for the case A//j > there exists a threshold value of the 
chemical potential (for fixed coupling constant Gs), or of the coupling con- 
stant Gg (for fixed chemical potential), below which we have only the trivial 
solution A = 0, and above which there is one unique minimum at A > 0. 
The transition to the color broken phase is therefore of second order in this 
case. On the other hand, for km = 0, the function /(A^) has a logarithmic 
singularity for A^ — > 0, and therefore the nontrivial solution exists for all 
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values of the chemical potential and of the coupling constant. In equations, 
if Ajji — the behavior for small A — > is expressed as follows: 

/(A^) = -^-^p^^flnA' + ... (4.48) 
A (X exp(- J , (4.49) 



where pp = ^l^q^ ~ M'^. In this case, Va has the Mexican hat shape for all 
values of /i* and Gs, and the system is always in the color symmetry broken 
phase. This behavior, which is very familiar from ordinary BCS theory [31], 
can be derived directly from QCD without using an effective quark theory 
[33]. This observation further supports our choice Air = in QM. 

In Fig. 8 we show the P — plots. The dotted line (r^ = 0) is very 
similar to the corresponding line (r^^ = 0) of Fig. 5, and shows that there 
is no phase transition from NM to normal QM. However, the softening of 
the QM EOS arising from diquark condensation leads to a first order phase 
transition, even for relatively small values of r^. For example, for rs — 0.1 the 
phase transition appears only at extremely high densities, but for rs — 0.15 
the phase transition sets in at 0.72 fm~^ in NM, and ends at 1.21 fm~^ in QM. 
The corresponding transition densities for = 0.2 are 0.57 fm~^ and 0.95 
fm~^. For the case = 0.25 the transition densities are already somewhat 
too low (0.42 fm~^ and 0.77 fm^^), and Vs = 0.3 is definitely too strong since 
QM would be the ground state for practically all densities. One should also 
note that, for any reasonable scenario, those parts of the QM curves which 
reflect the chiral phase transition are always below the NM curve. The chiral 
phase transition in QM is therefore not very relevant for the discussion of 
the NM QM phase transition, and one could set M = in QM from the 
beginning with practically identical conclusions. 
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In Fig. 9 we show the quark effective masses as functions of the density 
for the same values of as in Fig. 8, and in Fig. 10 we show the gaps as 
functions of the density. The behaviors of the curves M* (p) still reflect the 
chiral phase transition in QM, although the density region where M* varies 
rapidly is shifted upward compared to the case Vg — 0. Fig. 10 shows that 
the gaps in the present calculation are rather large. In the density region 
where QM becomes the ground state the gap is of the order of 150 MeV to 
250 MeV for those values of which give reasonable transition densities. 
This figure also supports the conclusion drawn from Fig. 8 that = 0.3 is 
already too strong, since the gap would exceed 300 MeV in the high density 
region, which is definitely too large. 

With regard to the dependence of the above results on the pairing strength, 
Tg, we would like to make two remarks. The first one concerns the regular- 
ization scheme. If we compare, for example, the present proper time scheme 
to the sharp 3-momentum cut-ofT scheme, a glance at Table 1 of Ref.[17] 
shows that the coupling constant Gt^ in the proper time scheme is almost 
three times as large as the one in the 3-momentum cut-off scheme. There- 
fore, for a given value of r^, the actual value of the coupling constant, Gg, 
in the proper time scheme is much larger than in the 3-momentum cut-off 
scheme, which partially explains why our gaps are larger than in previous 
works [7]. One also has to note that the most important contributions to 
the gap equation (4.45) come from the region \k\ ~ /x* - see Eqs.(4.46) and 
(4.47). In the high density region, this is of the same order as the UV cut-off 
in the 3-momentum scheme, and therefore the peak contributions to the gap 
equation are cut out artificially. On the other hand, the proper time scheme 
introduces a smooth cut-off function, and the gap can become quite large in 
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the high density region. 

The second remark is that the value = 0.51, obtained by fitting the 
nucleon mass in the pure quark-scalar diquark model, leads to unphysical 
results in QM. Here again we face the problem that the parameters which 
work for the nucleon and NM lead to unphysical results when used in QM. 
However, we beheve that for the case of the reason is more simple. The 
pure quark-scalar diquark model attributes the whole attraction in the 3- 
quark system to the scalar diquark correlations, whereas it is known that 
there are other important attractive effects, such as pion exchange [34, 35]. 
Indeed, it has been shown in Ref.[36] that the inclusion of pion exchange 
leads to a large reduction of in order to reproduce the experimental nucleon 
mass. The axial vector diquark channel also leads to some reduction of 
[9]. It therefore seems that = 0.51 is actually an overestimate, and if 
is determined by using a more complete model for the nucleon, it might also 
work in the QM calculation. 

Let us now illustrate the NM QM phase transition, taking the case 
Vg — 0.2 of Fig. 8 as an example. Since for a given // the actual ground 
state of the system is the state with the largest pressure, we see from Fig. 
8 (and the insert shown in Fig. 5) that below the NM saturation point the 
vacuum (VAC) is the ground state, and then the phase transitions to NM 
and QM occur successively. Since both phase transitions are of first order, 
the chemical potential is continuous, but the density is discontinuous across 
the transitions. In Fig. 11 we show the density of the ground state as a 
function of fi. The density jump for the NM— > QM transition is quite large, 
about 0.4 fm~^. Fig. 12 shows the constituent quark mass and the gap of 
the ground state as functions of This figure illustrates again that the NM 
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state shows only little tendency toward chiral symmetry restoration, and the 
transition leads to a QM state which is (almost) chirally symmetric, but 
where the color symmetry is strongly broken. In Figs. 13 and 14 we show 
the pressure and the energy density of the ground state as functions of the 
density. In the mixed phases, which are indicated by the dashed lines, the 
pressure is constant and the energy density increases linearly with density. 

5 Summary and conclusions 

In order to discuss the behavior of baryonic matter over a wide range of 
densities, one has to describe the single hadron, normal saturating nuclear 
matter, and high density quark matter at the same time. In this paper we 
have shown a consistent formulation and a feasible numerical treatment of 
these three aspects, which makes use of the quark-diquark picture for the 
single nucleon and the Hartree approximation for the many-body systems. 
We used the NJL model as an effective quark theory, since this model seems 
to be a unique candidate to provide a simultaneous covariant description of 
all three aspects of the problem. We paid special attention to the question 
of those conditions under which a phase transition from nuclear to quark 
matter becomes possible. 

For the description of normal saturating nuclear matter, there are two im- 
portant ingredients in the NJL model. First, one has to take into account one 
aspect of confinement physics, namely the absence of unphysical thresholds 
for the decay into quarks. In the NJL model this can be realized by intro- 
ducing an infrared cut-off in the framework of the proper time regularization 
scheme. Only if these unphysical thresholds are avoided, can one describe 
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normal saturating nuclear matter and stable hadrons in the medium. Sec- 
ond, the vector interaction leads to a repulsion between the nucleons which 
is inevitable for describing saturation. 

These two ingredients are very important in nuclear matter, but it is not 
clear a priori whether they simply can be taken over to quark matter. In 
this paper we presented arguments which indicate that it is preferable to 
use no infrared cut-off and no vector interaction in the high density quark 
phase. First, the use of an infrared cut-off would lead to hadronic poles 
for all densities in quark matter, and to a threshold behavior of the gap for 
small pairing strength and/or small chemical potential, which contradicts 
the results obtained on more general grounds. It is therefore preferable to 
use no infrared cut-off in quark matter. Second, the use of the same vector 
interaction as in nuclear matter would lead to uj meson poles in the high 
density region of quark matter. 

After the discussion of these points, we turned the problem of the phase 
transition from nuclear matter to normal quark matter. We found that there 
is no phase transition at all - i.e., nuclear matter is the ground state for all 
densities. We then took into account the effect of scalar diquark condensa- 
tion in quark matter, treating the 4-fermi coupling constant in this channel 
(the pairing strength) as a free parameter. We found that the quark matter 
equation of state is softened considerably, even for relatively small pairing 
strengths, leading to a first order phase transition from nuclear matter to 

^^We remind the reader that, concerning the equation of state of high density quark 
matter, the value of the infrared cut-off has no strong influence on the results, but the 
strength of the vector interaction is crucial. There would be no phase transition from 
nuclear matter to quark matter if one used the same strength of the vector interaction for 
the description of the two phases. 
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quark matter at a transition density which decreases with increasing pairing 
strength. That is, in our calculation there is a phase transition from nuclear 
matter to superconducting quark matter, but no phase transition to normal 
quark matter. The ground state of the system is then nuclear matter at 
normal densities, where chiral symmetry is strongly broken and color sym- 
metry is intact, and quark matter at high densities, where chiral symmetry 
is almost restored and color symmetry is strongly broken. We illustrated this 
picture by using an explicit numerical example in our model calculation. 

The results discussed in this paper for the equations of state of nuclear 
matter and quark matter show that, in many respects, these two systems 
have very different properties. In order to discuss the phase transition more 
quantitatively, it is important to account for these differences in the frame- 
work of an effective theory which is able to describe both phases consistently. 
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Appendices 



A Hadronization in the path integral formal- 
ism 

In this Appendix we present some formulae which are used in sect. 3.1. The 
Lagrangian density is given by (2.6), supplemented by (3.9). 

A.l Functional integration formulae 

To derive (3.10), (3.11), the integration over the quark fields is done by using 
the formula 



JVi/j Jvip exp (i + + ^^]) = exp 



(A.l) 



S S 

To derive (3.21), (3.22), we replace A ^tt^, A* ^tt^ in the in- 



diiJ* 



teraction part Sdi, and use the formula 



Jt^^J VA^ex.p (i [AD-^A* + A* J + J*A]) = exp ( 



S{iJ) 



iTrlnD"^ - J*DJ 



)■ 

(A. 



Finally, to derive (3.27), one integrates over the nucleon sources cp and 
by using a relation analogous to (A.l), i.e.. 



y" P^y" P^exp [i [0Giv0 + 0iV + iV0]) = exp ( 



-itT\nGN + NGnN 



(A.3) 



Note that the trace symbol Tr in (A.l) includes a trace in Nambu-Gorkov 
space, in contrast to the trace in (A.3), and therefore there is no factor | in 
the exponent of (A.3). 
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A.2 Derivation of Eq. (3. 18) 

The first term in (3.11) involves the propagator S, which is obtained by 
inverting (3.12). If we write 

then the diagonal elements are obtained as 
'-'11 — z : — : — '-'22 — 



while the non-diagonal elements are 

'S'12 = i Su'lb^aPaSo, S21 = i S22lb^*al3aSQ. 

The quantities Sq and 5*0 are given in momentum space by (3.17). When 
the above form for S is inserted into (3.11), the non-diagonal elements of 5* 
lead to terms which involve different numbers of A or A*, and these terms 
vanish after the integration over the diquark fields. Therefore the first term 
in (3.11) effectively becomes 

W/ A^iiA* \ ^ 

Then one uses the Nambu-Gorkov forms of $ and which are expressed 
similar to Eq.(2.1) in terms of and 0, and the relations CSqC~^ — Sq, 
CSqC"^ — So, where the symbol T (transpose) stands also w.r.t. the func- 
tional matrix structure, i.e., it implies x x'. In this way one arrives at 
Eq.(3.18). 
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A.3 Derivation of Eq. (3.23) 



Let us extract the term oc 00 from (3.21) and (3.22). The second exponential 
factor in (3.21) gives 



1 + ^Do(f)SQ(, 



exp ( — Tr In 

1 - (^o/^o) + ■ ■ ■ = 1 + i ^n^v (/> + 



(A.6) 
(A.7) 



where a color factor 3 comes from the color trace in (A.6), and Hn = iSoDo 
as in (3.25). Next, the factor in (3.22) gives to O 



exp {iSDi (A, A*)) exp {-iJ* D J)^ ^^^^^^ 
= 1 + i \Sdi (A, A*) exp (-i J* Dq J) 



J=J*=0 



+ ... 



(A.8) 



1 + ^0 



where A = 



+A5o^ {K^Pa) [K^a) ^So {Aa'Pa') (A^) ^qA* + . . . 

xexp {-iJ* Dq J)}j^j*^o + . . . , 
S S 



4> 



(A.9) 



A* 



A contraction is obtained if one A and one 



5{ij*y 6{ij)' 

A* hit the same factor {—iJ*DQj) = {iJ*){iDQ){iJ), which appears in the 
expansion of the exponent in (A.9). This contraction of A^ with A^' gives 
a factor iDgSa'a- Among the many contractions which emerge in (A.9), the 
ladder graphs appear as a subset besides self interaction graphs [13]. To 
obtain the ladder graphs from (A.9), one contracts A on the most left side 
with A to the nearest to the right of it, then the next A on the left side with 
A nearest to the right of it, etc. To write down the resulting series, we make 
the color structure more explicit by decomposing 

(A„,/3„,) (A:/5„)]^^ ^ c^tK^i = (-3 Pi?)"" + ^pS^^'^ A:A* , 
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where the color singlet and octet projection operators are given by: 



(0)a6 



kl 



ka^lh 



3(8)a6 
kl 



^kl^ab — -^Ska^lb 



Then the ladder-type contribution from (A. 9) is given by 



I — 

1 + -0 

3^ 



1 + i4>IlN Tn Hat^ 







(A.IO) 



where the quark-diquark T-matrix is given in (3.26). The sum of (A. 7) and 
(A.IO) is then 1 + i^Gjv^ = cxp {i(j)GN(t)^ + . . . , where the quark-diquark 
propagator in the color singlet channel, Gn, is given by (3.24). In this way 
one arrives at (3.23). 



B Explicit evaluation of the Dirac determi- 
nant 



The contribution of the quark loop to the effective action is (see Eq.(3.38)) 
— -TrlnS""^ = — / d^xVf, where the inverse quark propagator is a 2 x 2 
matrix in Nambu-Gorkov space and given by (3.36). By using the relation 

' A B 



Trln 



C D 



j =Trln(-SC + SDS-^A) 



and 01 = -diag (0, 1, 1) = -Ci, we obtain 



= ^ / ^Tr In (A^ d + (-^ -M + z/70) - M + z/70)) , (B.l) 

where we defined u = fi* to simplify the notation. For the evaluation of the 
Dirac determinant it is sufficient, because of rotational invariance, to consider 
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the vector p along the z-axis, and (B.l) becomes 

Ve^y ^Tr In (A^ C^ - pi + + + - 2M^7° + S^Pt't") , 

where in this Appendix we use the notation = and Ep = \Jp^ + M^. By 
using exphcit representations for the Dirac matrices, it is easy to calculate 
the Dirac determinant exphcitly. Because of isospin we get a factor 2, and 
one is left only with the color determinant: 



= 2z / -^In detc [(A^ C.-p^ ^1 + - 4EJ 



{2n) 



(B.2) 



By using Ci = diag (0, 1, 1), one can calculate the color determinant with the 
result 



Vf = 2i 



(2vr)' 



{2 In [{pI - el) {pI - e'_)] + In [{pi - el+) {pi - el_)] } , 



(B.3) 



where = {Ep ± z/)^ + A^ and el^ — {Ep ± i/)^. In this way one arrives at 
(3.38). 

It is convenient to separate the part which survives for A = from the 
rest. The latter contribution to the effective potential is given by the first 
line of (3.41). The former contribution becomes, after subtracting the zero 
density value, 

d'^p 

\4 



y.(A = 0)=6^/^ 



Vl - El 



-'Op PO ^Op 

where Eop — ^JMq + p^. This can further be spht into the pure vacuum loop 
which contributes to Vy^c given in (3.34), and the density dependent term 



, (B.4) 



Vq = Qt 



d'p [^^ pI - {Ep + ^ ^y^-{E,-uy 



(27r)^ L pI - El 



pI - El 
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= -6 



= -12 




J (27r) 



3 [{Ep + iy)-Ep + \Ep -iy\- Ep] 



(B.5) 



which is Eq.(3.40). 

We also note that the evaluation of the nucleonic contribution, (3.30), to 
the effective potential for NM proceeds similarly, by setting A = in (B.l) 
and replacing M — > Mjv, — > (Of course, there is no color factor 3 in 
this case.) Since the contribution for // = is subtracted, one is left with the 
density dependent piece (3.32). 

C The ou meson in nuclear and quark matter 

The expression for the the qq bubble graph in quark matter (Fermi momen- 
tum pp) for an external isoscalar vector field with q^ = (g°, 0) is given by 



We introduce a Feynman parameter (x) and make a shift k k — qx. Since 
q = 0, this shift affects only the time component, which is no problem 
since the proper time regularization scheme does not directly cut the loop 
momentum. The evaluation of the Dirac trace gives a term oc k'^k'^, for which 
we perform a partial integration, thereby picking up a surface term arising 
from the O function in (C.l). The result is no°(g) = and U'^{q) = g'm{q) 




(C.l) 



with 



U{q) = -48 g2 



/ 



d^k 




x{l — x) 



(27r) 



(A;2 + M2 - q^x{l - x)f 



(C.2) 
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1 



(C.3) 



/o k'^ + M'^- q^x{l - x) 

Here we performed a Wick rotation in ko. This expression is now regularized 
according to the prescription (4.42). (In order that the imaginary part is 
canceled for the case Air > 0, the finite surface term (C.3) has to be treated 
in the same way as the divergent part (C.2).) For the part (C.2) we introduce 

Euclidean polar coordinates (/c, 6) by \k\ = fcsinG, k^ = fccosO, where © 

Pf 

is restricted by < 6 < arc sin—. After performing the integration we 

k 

obtain for the piece (C.2), which we call H, 



X exp 



dr r 



t + M^ - q^x{l - x) 



(C.4) 



where 



G{y) = 1 - - I arc smy -ys^ L-y 



The T integration can also be performed, and the final result is 



n(g) 



X 




(\xxil — x) — 



m 1 + 



(C.5) 



where A ^ t + M"^ - qlx{l - x). For Air > 0, [. . .] in (C.5) behaves oc 
for A — > 0, and there is no imaginary part. The piece (C.3) is given by 



exp {-B/K 



SU{q)^^ Tdko f'dx 

TT'^ Jo Jo 



where B = k^ + + Pf — Qo ^(1 ~ 



exp 



B 
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D Numerical results for other choices of pa- 
rameters 

In this Appendix we present some numerical results for the QM EOS obtained 
by using values of Air or other than those used in the main text. 

First we consider the case where the same strength of the vector interac- 
tion as in NM {r^ = 0.37) is used also in QM, taking Air = in QM. The 
results are shown in Fig. 15. The dotted line here (r^ = in QM) is very 
similar to the dashed hne for r^^ = 0.37 in Fig. 5. (We note that in Fig. 5 
we used Air = 0.2 GeV.) It is clear from Fig. 15 that it is not possible to 
have a NM — > QM phase transition, even if one uses large values of r^. 

Next we consider the case where the same value of the IR cut-off as in 
NM (Air = 0.2 GeV) is used also in QM, taking = in QM. The results 
are shown in Fig. 16. The dotted hne here {vg = in QM) agrees with the 
dashed hne for r^, = in Fig. 5. We see that the quahtative behavior of 
the curves with increasing Vg is similar to the case Air = shown in Fig. 8, 
although the sensitivity to the pairing strength is somewhat stronger for 
the case Air = 0. (For given r^, the gaps are larger for the case of Air = 0.) 
The case of = 0.25, shown in Fig. 16, is actually very similar to the case 
Vg — 0.2 in Fig. 8, which was used in sect. 4.3 as an example to illustrate 
the NM — > QM phase transition. 
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Tables 





Nucleons and NM 


QM 


m [MeV] 

G^[GeV-2] 
Auv[MeV] 
Air [MeV] 

rs 


16.93 
19.60 
638.5 
200.0 
0.37 
0.51 


17.08 

19.76 
636.7 





free parameter 



Table 1: Parameters used for single nucleons, nuclear matter (left column) 
and quark matter (right column). The proper time regularization scheme is 
used in both cases. 
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Figure Captions 



1. Nucleoli mass as function of scalar potential for Air = 0.2 GeV and 
■^IR = 0. The parameter set for Air = 0.2 GeV is shown in Table 1, and 
that for Air = can be found in Table 1 of Ref.[17]. The quark-diquark 
threshold for the case Air = is also shown. 

2. Binding energy per nucleon as function of density for Air = 0.2 GeV 
and Air = 0. The parameter set for Air = 0.2 GeV is shown in Table 
1, while that for Air — can be found in Table 1 of Ref.[17]. 

3. The nucleon and quark effective masses in NM (solid hues), and the 
quark effective mass in normal (non-superconducting) QM (dashed 
line) are shown as functions of the density. The parameters used in 
the NM calculation are shown in the second row of Table 1, and the 
result for QM shown here refers to the same parameter set. (Note that 
the parameters r^, and rg are relevant only for the NM result shown in 
this figure.) 

4. Effective mass of the u; meson at rest (qr = 0) in NM as function of 
the density. The dashed hne shows the density dependence of 2M*. 
The parameters used in the calculation are shown in the second row of 
Table 1. 

5. P—fji plots for NM (sohd hne and insert) and normal (non-superconducting) 
QM. The two dashed hnes show the QM result for two values of r^, 
and the dotted lines show the massless quark gas EOS for two val- 
ues of the bag constant B. The parameters used in the NM calcula- 
tion are shown in the second row of Table 1, and the results for QM 

49 



shown here refer to the same parameter set, apart from the variation 
of the parameter r^^. (Both dashed hnes in this figure start at the point 
(//, P) = (1200MeV, 0), corresponding to zero density. For the case 
r^^ — 0.37 the pressure increases monotonously, while for r^^ = the 
pressure first increases slightly for very small densities, then decreases 
to negative values (unstable branch), and finally increases.) 

6. Effective mass of the u; meson at rest {q — 0) in QM, as function of 
the density, for the case Air = 0.2 GeV. The variation of the effective 
quark mass with density used in this calculation is shown by the dashed 
line in Fig. 3. The values of used in the calculation of the effective 
u;-meson mass are indicated for each line. The dashed line shows the 



7. Same as Fig. 6 for the case Air = 0. For small r^j, solutions exist only 
in the high density region. 

8. P — fx plots for NM (sohd hue) and QM for several values of the pairing 
strength in QM. The parameters used in the calculation are shown 
in the second and third rows of Table 1, and the values of in QM are 
indicated in the figure. A crossing of the NM and QM lines indicates 
a first order phase transition. The dotted line (r^ = 0) differs slightly 
from the result shown in Figure 5 for r^^ = 0, since for Fig. 5 the 
parameters shown in the second row of Table 1 were also used in the 
QM calculation. 

9. The effective quark masses in NM (sohd hue) and QM for several values 
of the pairing strength in QM as functions of the density. The cases 



density dependence of 2WM*2 + p|,. 
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1 to 6 indicated in this figure refer to different choices of Vg in QM as 
indicated in Fig. 8. The dotted fine (r^ — 0) differs from the result 
shown in Fig. 3, since for Fig. 3 the parameters shown in the second 
row of Table 1 were also used in the QM calculation. 

10. The gap in QM for several values of the pairing strength, , as a 
function of the density. The cases 1 to 6 indicated in this figure refer 
to the different choices of Vg shown in Fig. 8. 

11. The density of the ground state as function of the chemical poten- 
tial. The ground state of the system is either the vacuum (VAC), nu- 
clear matter (NM) or superconducting quark matter (QM). The pairing 
strength = 0.2 is used for QM. 

12. The quark effective masses (solid lines) and the gap (dashed line) of 
the ground state as functions of the chemical potential. The ground 
state of the system is either the vacuum (VAC), nuclear matter (NM) 
or superconducting quark matter (QM). The value of the gap is equal 
to zero for VAC and NM. The pairing strength Vg = 0.2 is used for QM. 

13. The pressure in the ground state as function of the density. The ground 
state of the system is either the vacuum (VAC), nuclear matter (NM) 
or superconducting quark matter (QM). (The vacuum corresponds to 
the point P = 0, p = 0.) The dashed lines indicate the mixed phases. 
The pairing strength = 0.2 is used for QM. 

14. The energy density of the ground state as function of the density. The 
ground state of the system is either the vacuum (VAC), nuclear mat- 
ter (NM) or superconducting quark matter (QM). (The vacuum corre- 
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spends to the point S = 0, p = 0.) The dashed hnes indicate the mixed 
phases. The pairing strength = 0.2 is used for QM. 

15. P — fi plots for NM (sohd hne) and QM for several values of the pairing 
strength, r^, in QM. The case shown here refers to Air = and r^^ = 
0.37 in QM. 

16. P — fi plots for NM (solid line) and QM for several values of the pairing 
strength, r^, in QM. The case shown here refers to Air = 0.2 GeV and 
To, = in QM. 
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